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Abstract 

We present a first-principles theory of the variation of magnetic anisotropy, K, with temperature, 
T, in metahic ferromagnets. It is based on relativistic electronic structure theory and calculation 
of magnetic torque. Thermally induced 'local moment' magnetic fluctuations are described within 
the relativistic generalisation of the 'disordered local moment' (R-DLM) theory from which the 
T dependence of the magnetisation, m, is found. We apply the theory to a uniaxial magnetic 
material with tetragonal crystal symmetry,Llo-ordered FePd, and find its uniaxial K consistent 
with a magnetic easy axis perpendicular to the Fe/Pd layers for all m and proportional to 
for a broad range of values of m. This is the same trend that we have previously found in LIq- 
ordered FePt and which agrees with experiment. This account, however, differs qualitatively from 
that extracted from a single ion anisotropy model. We also study the magnetically soft cubic 
magnet, the Fe^oPt^Q solid solution, and find that its small magnetic anisotropy constant Ki 
rapidly diminishes from 8 ^ueV to zero. K evolves from being proportional to m"^ at low T to 
near the Curie temperature. 

PACS numbers: 75.30.Gw, 75.10.Lp, 71.15.Rf,75.50.Bb,75.50.Ss 
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I. INTRODUCTION 



It is well-known that a description of magnetic anisotropy, K, can be provided once rel- 
ativistic effects such as the spin-orbit coupling on the electronic structure of materials are 
considered. Over recent years 'first-principles' theoretical work, based on relativistic density 
functional theory, has been quite successful in describing trends in K for a range of magnetic 
materials in bulk, film and nanostructured form^"^, e.g.^"^. These results can be fed into 
micromagnetic models of the magnetic properties to describe phenomena such as magneti- 
sation reversal processes in magnetic recording materials^^. There are also implications for 
electronic transport effects such as anisotropic magnetoresistance (AMR)^^. Until only very 
recently, however, the temperature dependence of K was assumed to follow that given by 
single ion anisotropy models developed by Callen and Callen and others over 40 years ago^^ 
This assumption has now been challenged by ab-initio electronic structure theory^^'^^. The 
unexpected dependence of the magnetic anisotropy of Llo-FePt, found in experiment^^"^^, 
to decrease in proportion with the square of the magnetisation, m(T'), is described well by 
the new theoretical treatments whereas the single ion magnetic anisotropy models fail. Ev- 
idently the itinerant nature of the electrons in metallic magnets like FePt is an important 
factor. 

In this paper we present a detailed description of our ab-initio theory for the temperature 
dependence of magnetic anisotropy. It involves a fully relativistic description of the electronic 
structure and hence includes spin-orbit coupling effects. The thermally excited magnetic 
fluctuations are accounted for with the, by now, well-tried, disordered local moment (DLM) 
picture. 

The study of temperature-dependent magnetic anisotropy has recently become particu- 
larly topical owing to extensive experimental studies of magnetic films and nanostructures 
and their technological potential. For example, fabrication of assemblies of increasingly 
smaller magnetic nanoparticles has great potential in the design of ultra-high density mag- 
netic data storage media.^^ If thermally driven demagnetisation and loss of data is to be 
avoided over a reasonable storage period, there is, however, a particle size limit to confront. 
A way of reducing this limit is to use materials with high magnetocrystalline anisotropy, K, 
since the superparamagnetic diameter of a magnetic particle is proportional to [ksT / K)^ , 
where /c^T is the thermal energy. Writing to media of very high K material can be achieved 
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by temporary heating. ^^'^^ K is reduced significantly during the magnetic write process and 
the information is locked in as the material cools. Modelling this process and improving the 
design of high density magnetic recording media therefore requires an understanding of how 
K varies with temperature. 

The temperature dependence of magnetic anisotropy in magnets where the magnetic 
moments are well-localised, e.g. rare-earth and oxide magnets, is described rather well 
by these single ion anisotropy models but it is questionable whether this will also be the 
case for itinerant ferromagnets.^^ Owing to its high uniaxial magnetocrystalline anisotropy 
(MCA) (4-10 10^ ergs/cm^ or up to 1.76 meV per FePt pair^^'^^) the chemically ordered LIq 
phase of cquiatomic FePt, has attracted much attention as a potential ultra-high magnetic 
recording density material. Indeed arrays of FePt nanoparticles with diameters as small as 
3 nm have been synthesiscd.^^'^^ For a uniaxial magnet like this, K is the difference between 
the free energies, and F^^'^'°^ of the system magnetised along (0, 0, 1) and (1, 0, 0) 

crystallographic directions. So for the first application of our theory we chose Llo-ordered 
FePt}^ Careful experimental studies of its fundamental magnetic properties. ^^"^^ find that 
over a large temperature range, K(T)/K{0) = (m(T)/m(0))" ,where n = 2 instead of n = 3 
as expected from the simple single ion anisotropy model. We found our ab-initio calculations 
to be in good agreement with this surprising result. Mryasov et al.^'^ independently examined 
the same issues with a different theoretical but complementary approach and drew the same 
conclusions. In this paper, after providing full details of the R-DLM theory of magnetic 
anisotropy, we explore whether this (m(T)/m(0))^ behavior is a general property of the 
MCA of Llo-ordered itinerant transition metal uniaxial magnets by investigating another 
important uniaxial magnetic material FePd. We also study the temperature dependence 
of a material which has cubic rather than the tetragonal crystal symmetry of Llo-ordered 
alloys, and which is magnetically softer, namely compositionally disordered FePt. 

In the next section we describe the temperature dependence of the magnetic anisotropy 
that emerges from classical spin models with single site anisotropy. We then review briefiy 
current approaches to calculating K from first-principles electronic theory of materials at 
T = OK. An outline of the 'disordered local moment' (DLM) picture of metallic magnetism 
at finite temperature precedes a description of its rclativistic generalisation. It is shown how 
the temperature dependence of the magnetisation, m{T), can be found. The key outcome 
from the R-DLM theory is the formalism for the magnetisation dependence of magnetic 
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anisotropy ab-initio. Applications to uniaxial Llo-FePd and cubic Fe^oPt^o follow and the 
final section provides a summary. 

II. SINGLE ION ANISOTROPY 

The MCA of a material can be conveniently expressed as K = K^g^{n) where the 
Kj^s are coefficients, n is the magnetisation direction and g^^s are polynomials (spherical 
harmonics) of the angles -3, (f, fixing the orientation of n with respect to the crystal axes, and 
belong to the fully symmetric representation of the crystal point group. As the temperature 
rises, K decreases rapidly. The key features of the results of the early theoretical work 
on this effect are revealed by classical spin models pertinent to magnets with localised 
magnetic moments. The anisotropic behavior of a set of localised 'spins' associated with 
ions sitting on crystalline sites, i, in the material, is given by a term in the hamiltonian 
Han — Si k-yg-y{ei) with Cj a unit vector denoting the spin direction on the site i. As the 
temperature is raised, the 'spins' sample the energy surface over a small angular range about 
the magnetisation direction and the anisotropy energy is given from the difference between 
averages taken for the magnetisation along the easy and hard directions. If the coefficients 

are assumed to be rather insensitive to temperature, the dominant thermal variation 
of K for a ferromagnet is given by K{T)/K{0) = {gi(e))T/ {giis))o The averages (■ ■ ■ )t are 
taken such that (e)T = i^iT), the magnetisation of the system at temperature T, and I is the 
order of the spherical harmonic describing the angular dependence of the local anisotropy i.e. 
I — 2 and 4 for uniaxial and cubic systems respectively. At low temperatures K{T)/K{0) ^ 
(m(T)/m(0))'(^+i)/2 and near the Curie temperature T^, K{T)/K{0) (m(T)/m(0))^ 

These results can be illustrated straightforwardly in a way which will be helpful for the 
development of our ab-initio theory. Consider a classical spin hamiltonian appropriate to a 
uniaxial ferromagnet. 

H ^ -^^JijCi-Cj - k^{no-eif (1) 

where Cj describes the orientation of a classical spin at site i and Jij and k are exchange 
and anisotropy parameters, uq is a unit vector along the magnetic easy axis. A mean field 

— * 

description of the system is given by reference to a hamiltonian h-Ci where the orientation 

— * — * 

of Weiss field h, i.e. h — hn, determines the direction of the magnetisation of the system and 
has direction cosines (sin 1? cos </?, sin 1? sine/?, cos'j?). Within this mean field approximation 
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the magnetisation m is 777.(7) = J eP(e)de where the probabihty of a spin being orientated 
along e is P(e) = e~f^^'"'^/Zo with Zq — J e~^'^"'^de. The free energy difference per site 
between the system magnetised along two directions rii and ri2 is 



K{T) = / ((no • e)'e-^'^"i-^ - {no ■ efe-f^^'^'-^)de (2) 



J 



If rii and ri2 are parallel and perpendicular to the magnetic easy axis Uq respectively then 



where g2 is the Legendrc polynomial (3(no ■ e')^ — l)/2. As a function of the magnetisation 
'm{T)/m{0), K{T)/K{0) varies quadratically near the Curie temperature Tc and cubically 
at low T. The same dependence can be shown for this simple spin model for the rate 
of variation of magnetic anisotropy with angle that the magnetisation makes with the 
system's easy axis, namely the magnetic torque^^ — —dK/d-d. 

III. AB-INITIO THEORY OF MAGNETIC ANISOTROPY 

Magnetocrystalline anisotropy is caused largely by spin-orbit coupling and receives an 
ab-initio description from the relativistic generalisation of spin density functional (SDF) 
theory.^ Apart from the work by Mryasov et al.^^ and ourselves^^, up to now calculations 
of the anisotropy constants K have been suited to T = OK only. Spin-orbit couphng effects 
are treated perturbatively or with a fully relativistic theory^'^^. Typically the total energy, 
or the single-electron contribution to it (if the force theorem is used), is calculated for two 
or more magnetisation directions, fii and 7x2 separately and then the MCA is obtained from 
the difference, A£J. A£J is typically small ranging from meV to /xeV and high precision 
in calculating the energies is required. For example, we have used this rationale with a 
fully relativistic theory to study the MCA of magnetically soft, compositionally disordered 
binary and ternary component alloys^^'^^ and the effect upon it of short-range^ and long 
range chemical order^^ in harder magnets such as CoPt and FePt. 

Experimentally, measurements of magnetocrystalline anisotropy constants of magnets can 
be obtained from torque magnetometry^^. From similar considerations of magnetic torque, 
ab-initio calculations of MCA can be made. There are obvious advantages in that the MCA 
can be obtained from a single calculation and reliance is not placed on the accurate extraction 




(3) 
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of a small difference between two energies. In particular the torque method has been used 
to good effect by Freeman and co-workers^^ in conjunction with their state-tracking method 
to study the MCA of a range of uniaxial magnets including layered systems. 

If the free energy of a material magnetised along a direction specified by n = ( sin cos 93, 
sin'i?sin(^, cosi? ) is F^^\ then the torque is 

r<»)^-^. (4) 

The contribution to the torque from the anisotropic part of F^") leads to a direct link between 
the gap in the spin wave spectrum and the MCA by the solution of the equation^" 

f = 7(nAfr")). (5) 

where 7 is the gyromagnetic ratio. Closely related to T^") is the variation of F^") with 
respect to and (p, i.e. T^{^,ip) = — and T^{-)!},ip) = — ^f^- As shown by Wang et 
al.^^, for most uniaxial magnets, which are well approximated by a free energy of the form 

= Fiso + K2 sin' ^ + K^ sin^ 1?, (6) 

(where K2 and K4 and magnetocrystalline anisotropy constants and Fiso is the isotropic 
part of the free energy), T^(i? = 7r/4, 9? = 0) = -{K2 + K^). This is equal to the MCA, 
AF = F^^'^'^^ — F^^'^'^\ For a magnet with cubic symmetry so that 

F^^^ ^ Fiso + Ki (sin^ i? sin' 2^p + sin' 2^) , (7) 

a calculation of T^{'d — 7r/2, cp — n/S) gives —Ki/2. In this work we present our formahsm 
for the direct calculation of the torque quantities T^{'&,ip) and T^i^-djip), and hence the 
MCA, in which the effects of thermally induced magnetic fluctuations are included so that 
the temperature dependence is captured. 

In our formalism the motion of the electron is described with spin-polarised, relativistic 
multiple scattering theory. An adaptive mesh algorithm^^ for Brillouin zone integrations 
is used in the calculations to ensure adequate numerical precision for the MCA to within 
0.1 //eV^''^. Since we characterise the thermally induced magnetic fluctuations in terms 
of disordered local moments, we now go on to describe this picture of finite temperature 
magnetism. 
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IV. METALLIC MAGNETISM AT FINITE TEMPERATURES - 'DISORDERED 
LOCAL MOMENTS 

In a metallic ferromagnet at T = OK the electronic band structure is spin-polarised. 
With increasing temperature, spin fluctuations are induced which eventually destroy the 
long-range magnetic order and hence the overall spin polarization of the system's electronic 
structure. These collective electron modes interact as the temperature is raised and are 
dependent upon and affect the underlying electronic structure. For many materials the 
magnetic excitations can be modelled by associating local spin-polarisation axes with all 
lattice sites and the orientations {e} vary very slowly on the time-scale of the electronic 
motions. These 'local moment' degrees of freedom produce local magnetic flelds on the 
lattice sites which affect the electronic motions and are self-consistently maintained by them. 
By taking appropriate ensemble averages over the orientational configurations, the system's 
magnetic properties can be determined. 

This 'disordered local moment' (DLM) picture has been implemented within a multiple- 
scattering (Korringa-Kohn-Rostoker, KKR)^^"^^ formalism using the first-principles ap- 
proach to the problem of itinerant electron magnetism. At no stage does it map the many- 
electron problem onto an effective Heisenberg model, and yet it deals, quantitatively, with 
both the ground state and the demise of magnetic long-range order at the Curie temperature 
in a material-specific, parameter-free manner. It has been used to describe the experimen- 
tally observed local exchange splitting and magnetic short-range order in both ultra-thin Fe 
films^^ and bulk Fe, the damped RKKY-hke magnetic interactions in the compositionally 
disordered CuMn 'spin-glass' alloys^^ and the quantitative description of the onset of mag- 
netic order in a range of alloys^^'^^. In combination with the local self-interaction correction 
(L-SIC)^^ for strong electron correlation effects, it also gives a revealing account of magnetic 
ordering in rare earths'^'^. Other applications of the DLM picture include dilute magnetic 
semiconductors^^ and actinides^^. Short-range order of the local moments can be explicitly 
included by making use of the recently developed KKR- nonlocal-CPA (KKR-NLCPA)^^'^! 

We now briefly recap on how this 'disordered local moment' (DLM) picture is imple- 
mented using the KKR-CPA and how a ferromagnetic metal both above and below the 
Curie temperature can be described. Our main objective in this paper is to explain its 
relativistic extension and show how this leads to an ab-initio theory of the temperature 
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dependence of magnetic anisotropy when relativistic effects are explicitly included. 



V. RELATIVISTIC DISORDERED LOCAL MOMENT THEORY 
A. General framework 

The non-relativistic version of the DLM theory has been discussed in detail by Gyorffy 
et al.^^'^'^ Here we summarise the general framework and concentrate on those aspects which 
are necessary for a description of magnetic anisotropy. The starting point is the specifi- 
cation of Q^'^\{e}), the 'generalised' electronic grand potential taken from the relativistic 
extension of spin density functional theory (SDFT)^'^^. It specifies an itinerant electron 
system constrained such that the local spin polarisation axes are configured according to 
{e} = {ei, 62, . . . , c/v} where N is the number of sites (moments) in the system. For magnetic 
anisotropy to be described, relativistic effects such as spin-orbit coupling upon the motion 
of the electrons must be included. This means that orientations of the local moments with 
respect to a specified direction h within the material are relevant. The role of a (classical) 
local moment hamiltonian, albeit a highly complicated one, is played by Q^''\{e}). Note 
that in the following we do not prejudge the physics by trying to extract an effective 'spin' 
model from O*^"^ ({e}) such as a classical Heisenberg model with a single site anisotropy term. 

Consider a ferromagnetic metal magnetised along a direction n at a temperature T where 
the orientational probability distribution is denoted by P^") ({e}) , and its average 

(e^) = J... J eiP(") ({e}) dei . . . dejv = n . (8) 

is aligned with the magnetisation direction n. The canonical partition function and the 

probability function are defined as 

Z(^) = J- J e-^"'''«^» dei . . . rfejv , (9) 

and 

P'-' ({e}) - , (10) 

respectively. The thermodynamic free-energy which includes the entropy associated with 
the orientational fluctuations as well as creation of electron-hole pairs, is given by 

F^") = -ilnZ(") . (11) 
8 



By choosing a trial Hamiltonian function, Q^"^ ({e}) with Z^"^ 
/.../e-^^o''({e})rfe^...dejv, 



g-/3nf)({e}) 





and Fq"-* = — -^InZQ"^ the Feynman-Peierls Inequality^^ imphes an upper bound for the free 
energy, i.e., 

where the average refers to the probabihty Pq"^ ({e}). By expanding Qq^^ ({e}) as 

({e}) = E ^'^'^ (^^) + I E (^^' + ■ ■ ■ ' (14) 
the 'best' trial system is found to satisfy 

/ ,^,\ 



pr m) = ^--^ (12) 

^0 



(ni»))°_-(sj<»))° = (n!,»')_ -(Jif) , (15) 
(n<«i>° - - = (iff . - (jif )° , (16) 

and so on, where ( )g. or ( )g. denote restricted statistical averages with ej or both Cj 
and kept fixed, respectively. For example, for a given physical quantity, X, 

^ /. . ./ X(^) ({e}) P(") ({e}) ciei . . . de,_ide,+i . . . de^ 



(17) 



with 



(eO = y. . .y Po^") ({e}) ■ ■ ■ ■ ■ ■ (^ejv . (18) 

The relationship, 

(X(")> = j {X^%^ pP (e,) de, , (19) 

is then obviously satisfied. 
B. Mean-field theory 

In this case we take a trial system which is comprised of the first term in Eq. (14) only, 

({e}) = E^f^^O ■ (20) 
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Therefore, the quantities defined above reduce to 

Zf^ He-^'^-'^^^') de,...den ^ Y[zf , = j ^^^^^ de, , (21) 

i i 

pt m) = n^f ^ (^^) ' ^ (^^) = '^^^ ' (22) 

i 

and 

Ff^--\Y.^^^f^ ■ (23) 
In order to employ condition (15) the following averages have to evaluated, 

^?)' = E n / (^0 (^.) E / ^ (^0 Pf" (e.) de, , (24) 

) + E l4^{e,)PP{e,)de,, (25) 



ei 

consequently, - hf^V = ouf^ (e^ - /a;f ^ (e'J i^^"^ (?,) cie', and Eq. (15) implies. 



(fi 



;f (e,) - (nr-)>! = I (uP (e;) - (0(^)>°J (e;) de; = const. , (26) 
which leads to the relationship 

(e.) = (n(^)>° . (27) 

The free energy is now given by 

Fin) ^ ^^(n) >% 1 ^ y /^(") (ei) In (e,) de, . (28) 

This is the key expression for our subsequent development of the magnetic anisotropy energy. 
(In the following we shall omit the superscript from the averages.) 

C. The Weiss field 

To proceed further one can expand u^^^ (e) in terms of spherical harmonics, 

(e) = ^^;%r,^(e) , (29) 

£,m 
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where the constant term, -^^oJ^i^QQ, does not enter the statistical averages and therefore can 
be taken to be zero. The coefficients, (jjf'^^ (cj), can obviously be expressed as 



(30) 



Keeping the leading term only, I — 1, with Yi,±i (cj) = \/ ^ {&x =t ^Cy) ) ^i,±i (ci) 



x-e^, a;^-"^ (cj) can be written as ujf (cj) = /i^-"'' • Cj with 



Furthermore, 



Zf^ = y" exp {-^hf ■ e,) de, 



(31) 



sinh/3/i|"^ , 



where /i,- 



(n) 



(n) 



, and the probability distribution is 



(32) 
(33) 



Pf^ (e,) 



i5K 



(n) 



An sinh ph. 



(n) 



exp ( —(3h\^'' ■ e 



(34) 



Thus the average alignment of the local moments, proportional to the magnetisation, is 



-> in) 1^ \ 



I3K 



(n) 



47rsinh j3h 



(n) 



Ci exp 



(35) 



from which rhf^^ — m\"'^ h^f^ and 



m. 



in) 



d\nZ 



(n) 



1 



(n) ' 



(-/JAf) 



(36) 



follow, where L{x) is the Langevin function. Since in the ferromagnetic state, fhi 
we finally can write the Weiss field, ^-"^ = h\^^ n as 



n 



(37) 



Note that an identical Weiss field h^"^^ associated with every site corresponds to a de- 
scription of a ferromagnetic system magnetised along n with no reference to an external 
field. 
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D. Averaging With The Coherent Potential Approximation 




Gyorffy et al.^^ we follow the strategy of the Coherent Potential Approximation (CPA)^^ as 
combined with the KKR method^^. The electronic charge density and also the magnetisation 
density, which sets the magnitudes, {/x}, of the local moments, are determined from a self- 
consistent field (SCF)-KKR-CPA'^^ calculation. For the systems we discuss in this paper, 
the magnitudes of the local moments are rather insensitive to the orientations of the local 
moments surrounding them^°. We return to this point later. 

For a given set of (self-consistent) potentials, electronic charge and local moment mag- 
nitudes {/li}, the orientations of the local moments are accounted for by the similarity 
transformation of the single-site t-matrices^^. 



where for a given energy (not labelled explicitly) (z) stands for the ^-matrix with effective 
field pointing along the local z axis^^ and R (cj) is a unitary representation of the O (3) 
transformation that rotates the z axis along Cj. In this work (z) is found by considering the 
relativistic, spin-polarised scattering of an electron from a central potential with a magnetic 
field defining the z-axis^^. Thus spin-orbit coupling effects are naturally included. 

The single-site CPA determines an effective medium through which the motion of an 
electron mimics the motion of an electron on the average. In a system magnetised along a 
direction n, the medium is specified by t-matrices, t\"^ , which satisfy the condition 



where the site-diagonal matrices of the multiple scattering path operator^ are defined as. 



U{ei)^R{ei) t,{z)R{ei) 



+ 



(38) 




(39) 




(40) 



-1 



(41) 



and 




(42) 
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In the above equation, double underlines denote matrices in site and angular momentum 
space, t^"^ is diagonal with respect to site indices, while stands for the matrix of structure 
constants^^'^^. Eq. (39) can be rewritten by introducing the excess scattering matrices, 

-1 



Xf' (e.) 



(n) 



-in) 



(43) 



in the form 



Xf) (^A ^ = . 



5(n) 



Thus, for a given set of Weiss fields, h[^\ and corresponding probabilities. 



>(«) 



in) 



4:77 sinh ph. 



(44) 



(45) 



Eq. (44) can be solved by iterating together with Eqs. (43) and (42) to obtain the matrices, 
i^^}. The integral in Eq. (44) can be discrctizcd to yield a multi-component CPA equation 
which can be solved by the method proposed by Ginatempo and Staunton'^'^. Care has to 
be taken, in particular for low temperatures where P^^^ (cj) is a sharply structured function, 
to include a large number and/or an adaptive sampling of the grid points. 



E. Calculation of the Weiss field 



In the spirit of the magnetic force theorem we shall consider only the single-particle 
energy (band energy) part of the SDFT Grand Potential as an effective 'local moment' 
Hamiltonian in Eq. (37), 

n(^) ({e}) -- j de fpD {e- i/(^)) 7V(-) {e- {e}) , (46) 

where is the chemical potential, fpo {^'-i ^^^'') is the Fermi-Dirac distribution, and 
A^*^"-* {e; {e\) denotes the integrated density of states for the orientational configuration, 
{e}. For a 'good local moment' system such as many iron and cobalt alloys, this frozen 
potential approximation is well-justified and discussed shortly. 

The Lloyd formula^^ provides an explicit expression for A'"^") (s; {e}) , 

A^(^) {e- {e}) = ATo {e) - ^ Imlndet (f^ {e- {e})"^ - {e)) , (47) 
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with A^o (s) being the integrated DOS of the free particles. The integrated DOS can be 
further decomposed into two more pieces 

7V(-) {e; {e}) = No (e) + TVf ) (e) + ATV^") {e; {e}) , (48) 

where 

ATf ) (£) = i Im In det r^") (£) = -- Im In det (t^^^ {e)~^ - G (e)) , (49) 

7j- — c 7j- \-c — / 

is independent of the configuration, and can be written as an integral over reciprocal wave- 
vector k-space 

ATf ) (£) = i Im y In det {i^^ (e)"' - Gq (k, s)) dk , (50) 

while 

AiV(") (e; {e}) = ~ Imlndet - ^ (e)"' - |(£; {e})"') r^") (e)) , (51) 

is the only configuration dependent part of N^'^^ {e; {e}) . Decomposing (e) into a site- 
diagonal, r^(") (s), and a purely site-off-diagonal term, t°J-"'^ (s), 

rf") (e) = {rg is) S,j} , (£) = {rgj. (6) (1 - S,,)} , (52) 

A?"*^'*) (e; {e}) can further be evaluated as 

A7V(") (£;{e}) = --ImlndetM^") (e; {e})-- Imlndet f / - X^^^ (£;{e})r°(") (e)) , (53) 

TT TT \— — — ^ / 

where 

M(")(e;{e}) =1+ (|(^;{e})-^-|f (e)-^)if )(e) (54) 
^ {Mf\s;e,) S,j} , (55) 

with the matrices, Mf\e;ei) = (£, ei)"\dcfincd in Eq. (41). K^^H^'Ae}) = 

^2^"'^ {s;ei) where X-"^ (£, ej), is defined in Eq. (43). Therefore, 

AAr(^) (e; {e}) = -i V Im In det Mf ^ (e; - ^ Im In det f / - X^") {e; {e}) t°^^^ {e)] . 

(56) 

Since, Mj-'^-' (e;ei) depends by definition only on the orientation ej, the restricted average, 
( of the first term simply equals to — Mm In det M^"^ (e\ e») , while the second term 
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requires more care. Namely, 

Tr In (^l - ^(") (£; {e}) z^^^^ (e)) = J] ^Tr (^^(") (e; {e}) z^^^^ (e)) (57) 

fc=i 

= E ^ E tr(2f!fU.;?J42.Js)Xlf(.;e.J...4«'(.;SJ4»)J.)) . 

(58) 

with tr( ) denoting a trace in angular momentum space only. Prom the single-site CPA con- 
dition, Eq. (44), it follows that the restricted average of the term k = 2 identically vanishes 
and that for higher order k terms the only elements in the sums which contribute are those for 
which each of the indices ii {I — 1, ... ,k) occurs at least twice. These backscattering terms 
are neglected in the single-site CPA averages. A useful property of the configurationally 
averaged integrated density of states given by the CPA 

(Ar(") (e; {e})> = (e) - \ Imlndet (ff (s)"^ - (s)) 

- - Im ^ ^In det Mf ^ (e; e^)^ (59) 

i 

is that it is stationary with respect to changes in the t-matrices, t^") (e), which determine 
the effective CPA medium. Indeed this stationarity condition can be shown to be another 
way of expressing the CPA condition^^. We will use this shortly in our derivation of a robust 
expression for the calculation of the MCA. 

The partially averaged electronic Grand Potential is given by 

(Q(")>. =-y" &/Fo(£;i/("))A^f)(£) + i y" ImlndetMf (£;eO , 

+ ^ ^ y"d£ (e; z/(")) Im (|lndet Mf (e; e,)^ , (60) 



and the Weiss field, h\ , can be expressed, using Eq.(37), as 

hf = Ay (e, -n) y dejFD (e; ^ Imlndet [e-e, 



dei . (61) 



The solution of Eqs.(61) and (35) produces the variation of the magnetisation m|"'' with 
temperature T with m^^"^ going to zero at T = T^'^\ When relativistic effects are included, the 
magnetisation direction n for which Ti"^ is highest indicates the easy direction for the onset 
of magnetic order. We can define a temperature range ATaniso — T^^"^ — T^^^"^ where rig and 
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rih are the system's high temperature easy and hard directions respectively, which is related 
to the magnetic anisotropy of the system at lower temperatures. Indeed an adaptation of 
this approach to systems such as thin films in combination with T — OK calculations may 
be useful in understanding temperature-induced spin reorientation transitions.^^ 



VI. THEORETICAL FORMALISM FOR THE MAGNETISATION DEPEN- 
DENCE OF MAGNETIC ANISOTROPY AB-INITIO 

In the ferromagnetic state, at temperatures more than ATaniso below the Curie temper- 
ature, the magnetic anisotropy is given by the difference between the free energies, F^'^\ for 
two different magnetisation directions, ni, n2, but the same magnetisation m and therefore 
the same values of the products of the Weiss field magnitudes with 13, i.e. (3h[^^^ = I3h[^^\ 
Within our DLM theory this means that the single site entropy terms in Eq.(28) for each 
magnetisation direction cancel when the difference is taken and the magnetic anisotropy 
energy MCA can be written 

AF(ni, na) = {fl^^'^) - {fl^^'^) . (62) 

This becomes AF(ni,n2) ~ - / de fpo {e;u^^^'^) ((M"i) (e)) - {N^^^'> (e))) where a small 
approximation is made through the use of the one chemical potential Using Eqs.(59) 

and (50) this can be written explicitly as 

AF(ni, na) = - y de fpo {e; i/^^^^) j Indet (l + t^^'^ {e)~^ - t^^^ {e)~^ T^'{e, k)) dk 
-Y^j defpn J [i^^^^^ (e,) Indet Mf^ e,) - (e,)lndetMf) (£;e,)] de, 

i 

(63) 

As well as ensuring that the Brillouin zone integration in the above equation is accom- 
plished with high numerical precision^'^^, care must also be taken to establish accurately 
the two CPA media describing the system magnetised along the two directions, rii and ri2. 
Eq. (44) has to be solved to high precision in each case. These steps were successfully taken 

and tested for our first application on LIq FePt^^. We have found however that a less 
computationally demanding scheme for extracting the temperature dependence of the mag- 
netic anisotropy can be derived by consideration of the magnetic torque. It is sufficiently 
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robust numerically to be applicable to a range of magnetic materials, whether hard or soft 
magnetically and in bulk, film and nanoparticulate form. 



A. A TORQUE-BASED FORMULA FOR THE MAGNETISATION DEPEN- 
DENCE OF MAGNETIC ANISOTROPY 



We return again to Eq.(28) for the expression for the free energy F*^") of a system mag- 
netised along a direction n — ( sin i? cos </?, sin i? sin (f, cos -d ) and consider how it varies with 
change in magnetisation angles -d and </?, i.e. — — ^f^, — — ^g^"^ . Since the single 
site entropy term in Eq.(28) is invariant with respect to the angular variations we can write 



d 



>(") 



(64) 



By using Eqs.(59) and (60), together with the stationarity of the CPA integrated density of 
states to variations of the CPA effective medium, wc can write directly 

dPP (e,) 



E 



IndetMf^ {e]ei)dei 



(65) 



According to the form of P/"'* (e^) given in Eq. (45) the principal expression for the magnetic 
torque at finite temperature is thus 



Ci 1 IndetM^ {e;ei)de, 



(66) 

For a uniaxial ferromagnct such as a LIq 3d-4d/5d transition metal magnet or a magnetic 
thin film, the performance of a single CPA calculation for appropriate values of the energy e 
is carried out at fixed values of the products (and therefore a chosen magnetisation m) 
and for the system magnetised along n — (sin7r/4, 0, cos7r/4). Subsequent evaluation of our 
torque expression, Eq.(66), i.e —T^{'& = 7r/4, (/? = 0) yields the sum of the first two magnetic 
anisotropy constants K2 and K4. Similarly —T^{'& = 7r/2,(p = tt/8) gives an estimate of the 
leading constant Ki/2 for a cubic system. In the appendix we derive T^{,p) for a magnet 
at T = OK and show that this is equivalent to Eq.(65) for the limit j3h —>■ 00, i.e. when 
T OK. 
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VII. THE CALCULATION OF THE TEMPERATURE DEPENDENCE OF THE 
MAGNETISATION, M{T), AND THE M(T) DEPENDENCE OF THE MAGNETIC 
ANISOTROPY 

In a first-principles implementation of the DLM picture, the averaging over the orien- 
tational configurations of the local moments is performed using techniques adopted from 
the theory of random metallic alloys. ^^''^^ Over the past 20 years, the paramagnetic state, 
onset of magnetic order and transition temperatures of many systems have been successfully 
described. All apphcations to date, apart from our earher study of FePt^^ and the cases 
described in this paper, have, however, neglected relativistic effects and have been devoted 
to the paramagnetic state where the symmetry turns the calculation into a binary alloy-type 
one with half the moments oriented along a direction and the rest antiparallel. Once rela- 
tivistic effects are included and/or the ferromagnetic state is considered, this simplicity is 
lost and, as is shown above, the continuous probability distribution, Pj^"'\ei), must be sam- 
pled for a fine mesh of angles and the averages with the probability distribution performed 
numerically. (Careful checks have to be made to ensure that the sampling of P/"^(ei) is 
sufficient - in our calculations up to 40,000 values are used.) Of course, in the paramagnetic 
state pj^"\(ii) = ^ so that a local moment on a site has an equal probability in pointing in 
any direction Cj. 

The local moments change their orientations, {ci}, on a time scale r long in comparison 
with the time taken for electrons to 'hop' from site to site. Meanwhile their magnitudes 
fluctuate rapidly on this fast electronic time scale which means that over times r, the 
magnetisation on a site is equal to /liCi. As a consequence of the itinerant nature of the 
electrons, the magnitude Hi depends on the orientations of the local moments on surrounding 
sites, i.e. /i, = fj,i{{ei}). In the DLM theory described above, /Xj = (/^i({e}))ei, so that the 
size of the local moment on a site is taken from electronic charge density spin-polarised along 
Cj and integrated over the site. An average is taken over the orientations {e} on surrounding 
sites and the local charge and magnetisation densities are calculated self-consistently from 
a generalised SDFT formalism and SCF-KKR-CPA techniques. 

Being a local mean field theory, the principal failure of the DLM to date is that it does not 
give an adequate description of local moment formation in Ni rich systems because it cannot 
allow for the effects of correlations among the orientations of the local moments over small 
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neighbourhoods of atomic sites. (In principle this shortcoming is now addressable using the 
newly developed SCF-KKR-NLCPA method^^'^^.) In this paper however we will focus en- 
tirely on 'good' local moment systems where the sizes of the moments are rather insensitive 
to their orientational environments. In these cases, for example, the self-consistently deter- 
mined local moments of the paramagnetic DLM state differ little from the magnetisation 
per site obtained for the ferromagnetic state. For example in paramagnetic DLM Llo-FePd, 
a local moment of 2.98 //s is set up on each Fe site whilst no moment forms on the Pd sites. 
For the same lattice spacings (c =0.381nm, c/a —1 - note we have neglected the deviation 
of c/a from ideal found experimentally) we find that, for the completely ferromagnetically 
ordered state of FePd at T = OK, the magnetisation per Fe site is 2.96/1^ and a small 
magnetisation of 0.32/iB is associated with the Pd sites. Likewise, the Fe^oPt^Q f.c.c. dis- 
ordered alloy has local moments of 2.92//^ on the Fe sites in the paramagnetic state whilst 
the ferromagnetic state has magnetisation of 2.93//^ and 0.22//^ on each Fe and Pt site 
respectively (a =0.385nm). We can therefore safely use the self-consistently generated effec- 
tive potentials and magnetic fields for the paramagnetic DLM state along with the charge 
and magnetisation densities for calculations for the ferromagnetic state below Tg. 
Our calculational method therefore is comprised of the following steps. 

1. Perform self-consistent scalar-relativistic DLM calculations for the paramagnetic state, 
T > Tc, to form effective potentials and magnetic fields from the local charge and mag- 
netisation densities (using typically the local spin density approximation, (LSDA)). 
This fixes the single-site t-matrices, t^^^ (z) 

2. For a given temperature and orientation, n = (sin -i? cosy?, sin -i? sin (/?, cos -i?) determine 
the hf^^^s ( and also the chemical potential i/*^"^ from Eq.(59)) self consistently: 

(a) for a set of Aj = Ph'f'^ determine the by solving the CPA condition, Eq. (44); 

(b) calculate new Weiss fields, Eq. (61); 

(c) repeat steps 2. (a) and (b) until convergence. [For a system where there is a 
single local moment per unit cell, this iterative procedure can be circumvented. 
A series of values of A(= /?/i*^"^) is picked to set the probabihties, P^'^\ei) (and 
magnetisations rh — (e), m — |m|). The Weiss field /i^^Hs then calculated from 
(61) and the ratio of /i" to A then uniquely determines the temperature T for 
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each of the initially chosen values of A and hence the temperature dependence of 
the magnetisation.] 

3. Calculate the torque, T^{cp) from Eq.(66) to give the magnetic anisotropy and also aver- 
age alignment of the local moments, rh^f^ (T), proportional to the total magnetisation, 
/iimj, from Eq. (35). 

4. Repeat steps 2. and 3. for a different direction, n' if necessary. 

(On a technical point: all integrals over energy ^ J de fpD i^*-"^) ■ ■ ■ ai'c carried out via a 
suitable contour in the complex energy plane and a summation over Matsubara frequencies^^. 
We use a simple box contour which encloses Matsubara frequencies up to f=:i 10 eV//i.) 

In the following examples for the uniaxial ferromagnets FePt and FePd, we have carried 
through steps 1. to 4. for i?i = j,ipi — where — —{K2 + K^, (T^ — 0) and also for 
^2 = f , = where = -^{K2 + §-^'4), {T^ = 0). For the cubic magnet disordered 
FePt we use = = as a numerical check, where both and T<^ should be zero, 
and 'i?2 = f ,<^2 = f where — and = 

VIII. UNIAXIAL MAGNETIC ANISOTROPY - FERROMAGNETS WITH 
TETRAGONAL CRYSTAL SYMMETRY 

Our case study in this paper is Llo-FePd and the trends we find are very similar to those 
found for Llo-FcPt. Wc carry out the steps 1-4 elaborated in the last section. Figure 1 shows 
the dependence of the magnetisation upon temperature. In this mean field approximation 
we find a Curie temperature of 1105K in reasonable agreement with the experimental value 
of 723K^^. (An Onsager cavity field technique could be used to improve this estimate, see 

without affecting the quahty of the following results for K.) Although the shortcomings 

3 

of the mean field approach do not produce the spinwave T2 behavior at low temperatures, 
the easy axis for the onset of magnetic order is deduced, n = (0, 0, 1) perpendicular to the 
layering of the Fe and Pd atoms, (not shown in the figure) and it corresponds to that found 
at lower temperatures both experimentally^^ and in all theoretical (T = OK) calculations^^. 
Figure 2 shows the magnetic anisotropy energy, AF((0, 0, 1), (1, 0, 0)) = —{K2 + K4) versus 
the square of the magnetisation. The same hnear relationship that we found for FePt^^ is 
evident, a clear consequence of the itinerant nature of the magnetism is this system. This 
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Temperature T(K) 

FIG. 1: The magnetisation of FePt versus temperature. The filled circles refer to a magnetisation 
along h = (1,0, 1). Tc is at 1105K with the easy axis, (0,0, 1). The full line shows the mean field 
approximation to a classical Heisenberg model for comparison. 

magnetisation dependence differs significantly from that produced by the single ion model, 
also shown in the figure. At T = OK, i^2 + -^4 is 0.335meV is in fair agreement with the value 
of 0.45 meV inferred from low temperature measurements on well ordered samples^^ (as with 
FePt, K decreases significantly if the degree of long-range chemical order is reduced). The 
value is also in line with values of 0.1 to 0.5 meV found by other ab-initio approaches 

From for both 'd = ^,ip = and ^ = ^,ip = the magnitudes of the MCA constants 
K2 and K4 are extracted and shown in Figure 3. The dominance of K2 is obvious but it is 
also clear that the dependence is followed closely by the total anisotropy, K2 + K4, and 
only approximately by the leading constant It is interesting to note that an anisotropic 
classical Heisenberg model leads to similar m dependence to K if treated within a mean field 
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FIG. 2: The magnetic anisotropy of FePd as a function of the square of magnetisation. The fihed 
circles show the calculations from the ab-initio theory, the full line KQ{m{T) / m{0))'^ and the dashed 
line the single-ion model function Kq < (72(e) >t / < 92(e) >o with Kq = -0.335meV. 

approach. To illustrate this point we show in Figure 4 the results of mean field calculations 
of K for a model with both single-ion and anisotropic nearest neighbor exchange,!. e. where 
the following hamiltonian is appropriate: 



(67) 



The full curve shows the single ion model results for the limit j" = J-*", which are also 
shown in Figure 3. At low T as m{T) 1, K(T)/K{0) has the familiar m'*^'"*"^)/^ form with 
/ = 1 for a uniaxial magnet. By introducing a small difference between and J-*", so that 
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FIG. 3: The magnetic anisotropy constants K2, K4 of FePd as a function of the square of magneti- 
sation. The filled circles and full line show the calculations from the ab-initio theory of the sum, 
the dashed line with filled diamonds describes K2 and the dotted line with squares shows 1^4. 

- Jll = O.OIJ^, K{T)/K{0) varies as m^. 

IX. CUBIC MAGNETIC ANISOTROPY - THE F.C.C. FesoPtso SOLID SOLU- 
TION 

Crystal structure is known to have a profound effect upon the magnetic anisotropy. Mag- 
netic anisotropy within a single ion anisotropy model decreases according to m'*^'+^)/^ at low 
T, (m ^ 1) and proportional to m' for small m at higher T. For materials with tetragonal 
symmetry, / = 2 as shown in Fig.2. On this basis a cubic magnet's K should possess an m 
dependence where / = 4, i.e. at low T and at higher T. In this section we show 
our results for the itinerant magnet, compositionally disordered Fe^oPt^Q. In this system 
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(M(T)/M(0)) squared 



FIG. 4: The magnetic anisotropy energy K calculated in a mean field approximation for a model 
of a uniaxial magnet which has both anisotropic exchange and single ion anisotropy. The full line 
shows results with single ion anisotropy only, k = 0.0027-*- and — jH = 0. The dashed line 
shows results for the same k and - jH = O.OIJ-^. 

the lattice sites of the f.c.c. lattice are occupied at random by either Fe or Pt atoms. The 
cubic symmetry causes this alloy to be magnetically very soft. Ordering into a tetragonal 
LIq structure of layers of predominantly Fe atoms stacked alternately with Pt layers along 
the (1, 0, 0) direction causes a significant increase of K. Okamoto et al^^ have measured K 
of FePt carefully as a function of compositional order and the trend, for T = OK, has been 
successfully reproduced in ab-initio calculations^^'^^. 

As with our earlier calculations for Llo-FePt^^ and FePd, this disordered alloy's mag- 
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netisation follows a similar T-dependence to that of a mean field treatment of a classical 
Heisenberg model. We find a Curie temperature of 1085K, again a mean field value which 
is in reasonable agreement with the experimental value of 750K^^. Figure 5 shows our cal- 
culations of the magnetisation dependence of the leading magnetic anisotropy constant Ki 
(Eq.7). At T = OK, AF((0,0, 1), (1, 1, 1)) = -Ki/3 is just 2.8/xeV (±0.1/ieV) , some three 
orders of magnitude smaller than the uniaxial MCA {K2 + K4) we find its Llo-ordered coun- 
terpart^^. Despite this small value we find that our method is robust enough to follow the 
magnetisation and T-dependence of Ki. Ki is determined from a calculation of where 
for 1? = 7r/2 and (p = tt/8 it equals Ki/2. As expected Ki decreases rapidly - Fig. 5 depicts 
Ki versus the fourth power of the magnetisation. At low T Ki varies approximately as rrJ 
whereas this dependence becomes m"^ for smaller M and higher T. Fig. 5 also shows the 
behavior of the single ion model for a cubic system for comparison. As with the uniaxial 
metaUic magnets already investigated, the ab-initio R-DLM results differ significantly. 



X. CONCLUSIONS 



We have shown that by including relativistic effects such as spin-orbit coupling into the 
Disordered Local Moment theory of finite temperature magnetism, the temperature de- 
pendence of magnetic anisotropy can be obtained. Magnetic anisotropy is determined via 
consideration of magnetic torque expressed within a multiple-scattering formalism. For uni- 
axial metallic magnets with tetragonal crystal symmetry, Llo-FePt and FePd, we find K 
to vary with the square of the overall magnetisation, m{T). This is at odds with what an 
analysis based on a single ion anisotropy model would find but in agreement with experimen- 
tal measurements for FePt. An interpretation in terms of an anisotropic Heisenberg model 
explains this behavior We suggest that this behavior is typical for high T^. transition 
metal alloys ordered into a tetragonal structure. We find the first anisotropy coefficient, K2 
to be dominant. We have also investigated the magnetic anisotropy of metallic magnets with 
cubic crystal symmetry which are very soft magnetically. In the example of the f.c.c. substi- 
tutional alloy, FesoPtso, the leading constant Ki decreases according to m" where n ranges 
between 7 and 4 as the temperature is increased. This behavior also difi^crs significantly 
from that of a simple single ion model. Application of this R-DLM theory of magnetism 
at finite temperature has been confined here to bulk crystalline systems. It also, however. 
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0.2 0.4 0.6 0.8 1 

Magnetisation to the fourth power 

FIG. 5: The magnetic anisotropy constant Ki of the cubic magnet Fe^oPt^Q as a function of the 
fourth power of the magnetisation, m^. The filled circles show the calculations from the ab-initio 
theory and the dashed line the single-ion model function < (72(e) >t / < 52(e) >o with = 

has particular relevance for thin film and nanostructured metallic magnets^^"^^ where it can 
be used to uncover temperature-induced reorientation transitions, e.g. Buruzs et al.^^ have 
recently applied the theory to Fe and Co monolayers on Cn(lll). Future possible applica- 
tions also include the study of the temperature dependence of magnetostriction, the design 
of high permeability materials and magnetotransport phenomena in spintronics. 
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APPENDIX A: TORQUE FOR T OK 

Concerning the MCA of a ferromagnet at T = OK, the relevant part of the total energy 

is 

F(") ^- J fpD {e; ly^^^) N^^^ (e) (68) 
where N^^^ (e) is the integrated density of states, 

Ar(") (£) = A^o (e) - ^ Im In det (t^ (n; e)'^ - {s)) , (69) 

and the inverse of the single site t-matrix is 



U {n; e)^R {n) [z] e) R {n)~ 



(70) 



Now R in) = exp iam{m-J_) where affi is the angle of rotation about an axis rh = {'zAn)/\'zAn\ 
and J is the total angular momentum. The torque quantity Ti^^ = —^§^, describing the 
variation of the total energy with respect to a rotation of the magnetisation about an axis 
u, is 

Tif-'l J defFD{e;iy^''^)lm^^[\ndet[mer'-g^{e))] (71) 
which can be written 

--\jd^f^D (e; ^("^) Im tr (^rf (e) ^ {R (n) t (n; e)-' R (n)+)) . (72) 



Since 



dRjn) 
dan 



i{J_-u)R (n) and 



dRjn)'' 



■ u)R{n), 



i 

For T^'ll), U ■ u) is just Jy(^,y 

Consider now our finite temperature torque expression, Eq.(65), i.e. 

T^ = ^^Jds fpn {s; i^^")) Yl J ^^^^ ^ ' (^4) 
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with a — or if and 



By definition, 



TT 



^ ' ^ In det M (£; e,) de, 

Ace 



lim 

Aa-»0 



Ids fpD {e- 5^ ^nn^ ^ | / i^^^'-""^ (e,) In det Mf> {e- e,) de, 



Pf) (e,) In det M^^'^^ (£;e,)c^e, 



Approaching T = 0, 



f*^"^ (ci) ^5{n-ei) , (e; 6^) ^ (e; n) , 
(cj) ^ 5 (n + An - Cj) , (s; e^) ^ (s; n + An) , 



,(n+An) 



while 



( 



In det Mf^ (£;ei) ^ Tr 

Therefore, 



TT 



which is equivalent to Eq.(72). 
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